Anomalous versus non-homogeneous diffusion in the ligand-binding equilibrium: 

when size matters 
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Measurements of protein motion in living cells and membranes consistently report transient 
anomalous subdiffusion which converges back to a Brownian motion with reduced diffusion con- 
stant at long times, after the subdiffusion regime. On the other hand, membranes are also non- 
homogeneous media in which Brownian motion may be locally slowed-down due to variations in 
lipid composition. Here, we investigate whether both situations lead to a similar behavior for the 
reversible ligand-binding reaction in 2d. We compare the (long-time) equilibrium properties ob- 
tained with transient anomalous diffusion to those obtained with slowed-down Brownian motion. 
We show that both processes increase the apparent affinity of the reaction. However, in the case of 
slowed-down Brownian motion, the affinity is maximal when the slowdown is restricted to a subre- 
gion of the available space. Hence, these two processes are different: size matters for slowed-down 
Brownian motion, not for transient anomalous subdiffusion. 
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Cell membranes are heterogeneous collections of con- 
tiguous spatial domains with various length and time 
scales (e.g. lipid rafts, caveolae) pQ, that spatially mod- 
ulate the lateral diffusion of proteins pi. This defines 
an non-homogeneous diffusion problem, with position- 
dependent diffusion constant [3 . On the other hand, nu- 
merous studies [1H6] have shown that the lateral diffusion 
of lipids, proteins and multimolecular assemblies in cells 
and membranes exhibit subdiffusive anomalous diffusion 
(whereby the mean squared displacement scales sublin- 
early with time, (r 2 (i)) cx t a with a < 1), a hallmark 
of diffusion obstruction by obstacles [7] and, more gener- 
ally, heavy-tailed jump-time distributions [3 [5]. How- 
ever, in the latter situation, the anomalous regime is 
usually only transient: at long times, the mean squared 
displacement crossovers back to normal (Brownian) dif- 
fusion, with a = 1 but a reduced apparent diffusion 
constant [6]. Such transient behaviors are for instance 
obtained when the density of obstructing obstacles is 
below the percolation threshold [5]. This asymptoti- 
cally slowed-down Brownian regime could be considered 
a mesoscopic (homogenized) representation of the under- 
lying microscopic anomalous subdiffusion. At first sight, 
both situations could thus be considered equivalent on 
long time scales. Regarding the biochemical reactions 
that take place in these media, the influence of both sit- 
uations has been relatively unexplored. Some results are 
available for transient reactions with transient anoma- 
lous subdiffusion - for the first binding times [TU] and for 
Michaelis-Menten enzyme reactions and substrate con- 
sumption |lf) . For position-dependent diffusion, the re- 
sults are even sparser [T2] . 

Here, we question the validity of the approximation 
that slowed-down Brownian motion could account for 
transient anomalous diffusion at long times. We study 



the equilibrium properties of a reversible reaction when 
diffusion is transiently anomalous due to obstacles (be- 
low the percolation threshold) or when normal space- 
dependent Brownian diffusion takes place. We consider 
the ubiquitous ligand-binding equilibrium 

L + R ^ C (1) 

feoff 

where L is the ligand, R its free receptor, and C the 
bound complex. Reversible reactions are expected to 
converge at long times to equilibrium, thus permitting 
the study of the influence of both anomalous diffusion 
and space-dependent diffusion not only on transient but 
also long-time equilibrium properties. Note that in clas- 
sical conditions, the equilibrium of Eq.Q is well known 
(assuming L ^> i?n, where Lrj an d Ro are the ini- 
tial concentrations in ligand and receptor, respectively): 
C oq = RoL /(Krj + L ). This defines so-called dose- 
response curves - the equilibrium amounts of C for 
increasing doses of ligand - with equilibrium constant 
Kd = fc on /fc ff , a measure of efficiency of the reaction. 

To simulate diffusion, we initially position L and R 
molecules (and possibly, obstacles) uniformly at random 
on a 800 x 800 2D square lattice with periodic bound- 
aries. Each lattice site is associated with a dif- 
fusion constant D{i,j) (all molecules here have identi- 
cal diffusion constants). D(i,j) = if contains 
an obstacle. At each time step St, every molecule is al- 
lowed to leave its current location with probability 
— 2dSt/(Sx) 2 D(i,j), where 5x is the lattice spac- 
ing and d the dimension. The destination site is cho- 
sen uniformly at random from the 4 nearest neighbors 
(i ±l,j± 1) and the molecule jumps to it, except if the 
destination site contains an obstacle, in which case the 
molecule remains in its initial location If the jump 
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of the molecule to its destination site results in the pres- 
ence of a (L,R) couple on the same lattice site, a binding 
event may occurs, i.e. the (L,R) couple is replaced by a 
single C molecule at the site, with probability p on . Fi- 
nally, at each time step, every C molecule can unbind, 
i.e. the C molecule is replaced by a (L,R) couple at the 
same site, with probability p a s- 
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FIG. 1. Simulations of reaction eq.Q in homogeneous con- 
ditions, D(i,j) = D, Vi,j. (A) Evolution of the bound frac- 
tion at equilibrium, C cq /i?o as a function of the relative lig- 
and concentration Lo/Kno, where Kdo is the value of Kd in 
the absence of obstacles and with reference diffusion constant 
Do- (p, D)=(0.0, D ) (black), (0.0, D /2) (orange) or (0.35, 
Do) (green). The inset shows the initial regime of the cor- 
responding mean-squared displacements (MSD). (B) Relative 
apparent equilibrium constant Kd/Kdo for D = Do and in- 
creasing p (Bl) or no obstacles but increasing reduction of the 
diffusion, 7=1 — D/Do (B2). The dashed line locates the di- 
agonal y = 1 — x. Other parameters: Ro = 100, St = 5x = 1, 
Pon = 0.1, Pog = 10~ 3 and D = 1, data are averages of 10 
runs. 

In a typical simulation, we start with R R molecules 
and Lq L, and run the simulation until the density of 
bound receptors C reaches a steady state, C cq . As shown 
in Figure [T]4. we obtain the typical dose-response curves 
for the bound fraction at equilibrium C eq /i?o- From this 
we can retrieve for all situations the apparent equilibrium 
constant Kd- The control behavior - no obstacles and 
spatially constant diffusion, D(i,j) — Do (black curve) 
- is not changed when the diffusion constant is lowered 
over the whole space to Do/ 2 (orange curve). Whereas 
the presence of obstacles strongly modifies the equilib- 
rium (green curve, obstacle density p — 0.35), yielding 
lower Kd- Fig[Tj3 displays Kd values for several obsta- 
cles density p (Bl) and several diffusion constant reduc- 



tions 7 in the absence of obstacles (B2). Clearly, for these 
parameters the regime is reaction-limited: one order of 
magnitude span for D does not influence Kd- Far from 
the percolation threshold (p — 0.41), Kd decays linearly 
with obstacle density as Kd/Kdo = 1 — P- This is a 
simple effect of the excluded volume occupied by the ob- 
stacles since it disappears if the concentrations are com- 
puted on the basis of the accessible space 1 — p, instead 
of the whole space (not shown). Therefore, hindered dif- 
fusion due to obstacles not only decreases the molecule 
mobility, it also increases the affinity (~ 1/Kd) of the 
reaction. 

The results presented so far are for homogeneous con- 
ditions, i.e. D(i,j) = D, But in living cells, the 
conditions are usually spatially non-homogeneous : spa- 
tial domains (lipid rafts, caveolae) give rise to position- 
dependent values of D; obstacle density can also be 
spatially-variable, as major obstacles can be clustered. 
In the following, we addressed this situation by restrict- 
ing the region of space where diffusion is modified to a 
central square patch of variable spatial extent. Diffusion 
is Brownian with coefficient D outside the patch and 
reduced within the central patch either by the presence 
of obstacles or by imposing a reduced diffusion constant 
D\ inside the patch. 

Before addressing the ligand-binding equilibrium 
eq.([T]) in these conditions, we first investigate diffusion 
effects alone. We used a central patch which surface area 
was 25% of the whole space and simulate the diffusion 
of diffusive non-reactive molecules until they reach equi- 
librium. Once equilibrium is reached, we perturb it by 
the addition of supplementary non-reactive molecules in 
the center of the patch and measure the characteristic 
time to reach a new equilibrium and the concentration 
of molecules inside the patch at this new equilibrium. 
Because reduced diffusion in the patch slows down the 
molecules, the characteristic time to converge back to 
equilibrium increases when diffusion is reduced in the 
patch, either because of obstacles (Fig. [2j\l) or directly 
by decreasing Di (Fig.[2j31). Both cases are remarkably 
similar in that respect. The main differences occur for 
the equilibrium concentrations. With obstacles, the equi- 
librium concentration within the patch slightly decreases 
with increasing obstacle densities (Fig. [2]A2) . The results 
obtained for space-dependent Brownian conditions how- 
ever strongly depart from the case with obstacles. We 
observe increasing accumulation of molecules within the 
patch when diffusion is decreased therein (Fig. [2j32) . We 
emphasize here these are equilibrium conditions. 

However surprising, this equilibrium effect can be di- 
rectly predicted in our system. Let us consider the Id 
case for simplicity, and a constant-by-part dependence 
of the diffusion constant D(x) = D 1: Vx £ [a,b] and 
D(x) — Do outside the patch [a, b]. Let us then con- 
sider a single molecule and let Tv(x,t) its probability to 
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FIG. 2. Diffusion in non-homogeneous conditions. Non- 
reactive molecules diffuse with coefficient Do (no obstacles) 
outside of a central patch that contains obstacles at density p 
(A 1-2) or in which the diffusion constant is reduced according 
to Di = Do("f — 1) (no obstacles) (Bl-2). The top row shows 
the characteristic time to reach equilibrium and the bottom 
one the molecule density at equilibrium in the patch. Data 
are normalized by the values obtained for p — in the patch 
(Al-2) or when 7 = (Bf-2). Full lines are guide to the 
eyes while the black dashed line in B2 shows the theoretical 
prediction eq. |7f . All other parameters are as in Fig. [T] 



be located at position x at time t [T5] : 

n(x, t + St) =q(x)ir(x, t) + n(x — Sx, t) (1 — q(x — 5x)) /2 



+n(x + Sx, t) (1 - q(x + Sx)) /2 



(2) 



where q(x) is the probability not to jump at each time 
step and is defined, using the jump probability j(x) — 
25t/(Sx) 2 D(x) above, asq(x) = l—j(x). Noting g(x, t) = 
(1 — q(x))w(x, t)/2 and developing g(x ± 5x,t) in series 
of x, one obtains at order 2 

7r(x, t + St)= q(x)ir(x, t) + 2g(x) + (Sx) 2 d xx g(x) 

= n(x, t) + (Sx) 2 d xx g{x, t) (3) 

Dividing by St and taking the limit St — > 0, one gets 

d t ir(x,t)=d xx {D(x)ir{x,t)) (4) 

where we used the expression of j(x) above to define 
D(x). Noting u(x, oo) the density of molecules at x at 
equilibrium, one expects from eq.Q 



D(x)u(x, oo) = H(D) 



(5) 



where T-l(D) is the spatial harmonic mean of the (space- 
dependent) diffusion function 



H(D) = 



D- 1 (x)dx 



Now, using the constant-by-part function for D{x) ex- 
pressed above, this yields u(x 7 oo) = 'H(D)/Di Va; G [a, b] 
and u(x,oo) — %{D)/Dq outside. The equilibrium con- 
centration inside the [a, b] patch thus equals that found 
outside the patch multiplied by D /Di. Hence the larger 
the slowdown of the Brownian motion inside the patch, 
the larger the accumulation inside it at equilibrium ex- 
plaining simulation results of Fig. [2j32. In the present 2d 
case, the total number of molecule in the patch Ni ns ^ e 
relates to total number -/V tota i, the surface fraction of the 
patch (j) of the total surface S and the diffusion constant 
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(6) 



(7) 

Eq. ^ yields a very good prediction of the increase of 
equilibrium concentration of Fig. [5J32. 

Such concentration changes are likely to modify the 
reaction in the non-homogeneous diffusion case. Using 
several sizes for the patch we computed the values of 
the apparent equilibrium constant Kjj. With anoma- 
lous subdiffusion due to obstacles (Fig. [3}\), the behav- 
ior reported in Fig. [T] is roughly conserved for all patch 
area fractions (j): Kr> decreases linearly with the obsta- 
cle density p far from the percolation threshold then in- 
creases back close to it. The amplitude of this decay in- 
creases with the patch area (as the total space occupied 
by obstacles increases). Considering the behavior ob- 
served in spatially homogeneous conditions (Fig. [l}31), 
one expects, far from the percolation threshold, Kd = 
K D o(l-4>)+KDoO--p)<f>, yielding K D /K D0 = l-pcj). In- 
deed, we found that the latter is a very good approxima- 
tion for the simulation results of Fig.[3]A. Close to the per- 
colation threshold, the increase of the first-collision time 
overcompensates by far the decrease of the re-collision 
times due to obstacles (not shown). As a result, the for- 
ward reaction rate fc on , and the apparent affinity, strongly 
decrease close to the threshold. 

The situation is quite different for space-dependent dif- 
fusion (Fig. [3j3) . Whatever the patch size, we also ob- 
serve that Kd decreases, but in this case, this is the 
result of the accumulation phenomenon reported above 
(Fig. [2^32) . Moreover, in this case, Kb exhibits a non- 
monotonous dependency with respect to the area frac- 
tion (/>, with a marked minimum. Hence, for a given 
value of diffusion reduction in the patch, 7 = 1 — Di/Dq, 
there exists an optimal value of the patch surface area 
that yields the highest affinity (Fig. [3j3, inset). Assum- 
ing that reactions in the patch and outside are both 
reaction-limited and separable, one can estimate the cx- 
tremum reached by C concentration using eq([7| for both 
R and L. From this, we derive a theoretical relation be- 
tween 7 and the optimal value of the patch surface, (p*: 
7 = 1 — 0*/(l — (j)*) (see Supplementary Material). The 
inset of Fig. [3b illustrates the good agreement with simu- 
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lation results. Therefore, non-homogeneous slowed-down 
Brownian systems exhibit a bonus to patchiness : the 
minimal value of Kr> is obtained when the patch occu- 
pies a subset of the available space. This is in strong con- 
trast with the equilibrium behavior obtained with tran- 
sient anomalous subdiffusion above, where the affinity 
increases monotonously with the patch area fraction. 
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FIG. 3. Simulations of reaction eq.|TJ in spatially non- 
homogeneous conditions. With transient anomalous subdif- 
fusion (A) or reduced Brownian motion (B) inside a central 
patch, the apparent equilibrium constant Kd/Kdo is shown 
as a function of the area fraction occupied by the patch (f> 
and the obstacle density p or the amount of diffusion reduc- 
tion in the patch, 7=1 — D 1 /Dq. The inset in (B) locates 
the value of <j> yielding maximal affinity (full circles) and the 
corresponding theoretical prediction (full line). All other pa- 
rameters are as in Fig. [I] 

In summary, we have shown that the equilibrium prop- 
erties of the ligand-binding reaction eq.([T]) differ when 
diffusion is slowed-down by anomalous subdiffusion due 
to obstacles or by a decrease of the diffusion constant due 
to, e.g., a change in lipid composition of the membrane. 
In both cases, slowing down the motion increases the ap- 



parent affinity of the reaction. But in the former case, 
the increase of affinity is due to the decay of the accessi- 
ble space, whereas in the latter, increased affinity is due 
to accumulation of molecules in the slowed-down zones. 
Moreover, in the case of anomalous subdiffusion, max- 
imal affinity is obtained when the obstacles are spread 
all over the available space, whereas in the case of a 
space-dependent diffusion constant, maximal affinity is 
reached when the slowed-down region occupies only a re- 
stricted area. It is very tempting to remark that, in living 
cell membranes, slowed-down regions (rafts) show a very 
patchy distribution, whereas bulky obstacles seem less 
systematically clustered in limited regions. Controlling 
the spatial extension of the areas with reduced lateral 
diffusion may thus be a way by which cells control the 
apparent affinity of the ubiquitous ligand-reaction bind- 
ing events. 
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